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This is a review of the parton cascade model (PCM) which provides a QCD-based 
description of nucleus-nucleus reactions at very high energy. The PCM describes 
the collision dynamics within the early and dense phase of the reaction in terms 
of the relativistic, probabilistic transport of perturbative excitations (partons) of 
the QCD vacuum, combined with the renormalization group flow of the parton 
virtuality. The current state of numerical implementations of the model, as well 
as its predictions for nuclear collisions at RHIC and LHC are discussed. 



1 The Parton Cascade Model 

The parton cascade modclE (PCM) was proposed in_L991 and developed further 
during the following years, especially, by GeigerHtjoB. Its aim is to describe 
the energy deposition, thermalization, and chemical equilibration of matter in 
ultra- high energy nuclear collisions, providing a full space-time picture of the 
event up to the moment when individual hadrons are formed. The model was, 
at least originally, not conceived as an "event generator" that would predict a 
full set of hadron momentum distributions in the final state. The parton cas- 
cade model code VNI, developed by Geiger, allows us to make such predictions 
because it contains the implementation of a hadronization scheme similar to 
those used in jet fragmentation algorithms. However, it should be clear that 
this runs somewhat counter to the original purpose of the PCM, namely, to 
explore the range of validity of perturbative QCD in nuclear reactions. Nev- 
ertheless, the predictions of the PCM with an added hadronization stage are 
experimentally useful. 

The conceptual basis of the PCM is the inside-outside cascade model 
of hadronic reactions, which implements the concept that new matter pro- 
duced in hadronic interactions at high energy is formed outside the intersecting 
world-tubes of the colliding hadrons. Bjorken's hydrodynamical model 13 was 
developed to describe the evolution of this newly formed matter in the central 
space-time region after thermal equilibration. Beginning in the mid-1980's it 
was realized that the deposition of energy into this region may be described 
in terms of ideas d^ived-Erom perturbative QCD (minijets) in the case of col- 
liding heavy nuclciH'Enl90. A computer code (HIJIKG.) incorporating some of 
these ideas was developed by Gyulassy and WangliMa 
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The partem cascade model combines these concepts into one unified scheme 
for the description of the space-time evolution of matter in nuclear reactions. 
Its three main ingredients are: 

1. The initial state is viewed as incoherent ensemble of partons determined 
by the nuclear parton distribution functions qf(x, Q 2 ) and g(x, Q 2 ), where 
the subscript / denotes t he quark flavor and g stands for the gluon distri- 
bution, x = p z / P is the longitudinal momentum fraction of the nucleon 
carried by the parton, and Q 2 is the parton "scale" or virtuality. Before 
any interaction occurs, Q 2 is generally taken as space-like. Our knowl- 
edge about the space-time structure of the nuclei before the collision and 
our limited information about the intrinsic transverse momenta of par- 
tons is then used to construct a model for the six-dimension phase space 
distributions of partons before the interaction: qf(r,p),g(r,p). The par- 
ton distributions are conveniently initialized at the scale Qq — (py) co ii of 
the average momentum scale of the primary parton-parton interactions. 

2. The time evolution of the parton phase distributions is governed by a rel- 
ativistic Boltzmann equation with "leading-log" improved lowest-order 
collision terms. Only binary interactions are allowed, but the final state 
can have (and generally has) more than two particles. As is well known, 
the higher-order improvement of the cross sections by means of the lead- 
ing logarithmic approximation is equivalent to the scale evolution of the 
parton distributions according to the DGLAP equation. Motivated by 
quantum mechanical considerations, the space-time picture of parton 
propagation before and after interactions is closely related to their off- 
shell propagation: the formation of a parton with virtuality scale Q takes 
a time T/(Q) rs U/Q. 

3. When the parton distributions become sufficiently dilute, they hadronize. 
In VNI the hadronization is described by a clustering algorithm, followed 
by the decay of excited hadrons. The transition is assumed to occur 
when the average virtuality of the partons falls below a critical value 
Q cr it ~ lGeV, because partons no longer scatter with sufficient energy. 

From a gradient expansion of the evolution equation for the parton Wigner 
distribution two equations can be derived Ej. The first equation 

P tl -^ jI F i (r, P ) = C i (r, P ) (1) 
describes the free propagation of partons which is intermittently modified by 
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collisions given by the binary collision terms C. The second equation 

p 2 -^F i (r,p) = S i (r,p) (2) 

describes the evolution of the parton distributions with respect to virtuality 
or "off-shellness" p 2 . This equation is a generalization of the usual mass- 
shell condition F(r,p) ~ 5(p 2 — m 2 ) to the case where the on-shell particle 
distribution cannot be defined. Si(r,p) describes the splitting of single off-shell 
partons into two partons of smaller virtuality. 

The two equations can be viewed as quantitative representation of Feyn- 
man diagrams of the type shown in Figure 1. The collision term d is rep- 
resented by the binary collision diagram contained in the box at the heart of 
the complex Feynman diagram, whereas the splitting term Si is represented 
by the branchings diagrams of the initial- and final-state partons.^ The Feyn- 
man diagram of Fig. 1 is finite only if the virtualities of all final-state partons 
are limited by some infrared cut-off /i 2 . In an isolated event /j 2 describes 
the hadronization scale, i.e. the virtuality scale below which partons can no 
longer be considered as approximately free, perturbation quanta. In a dense 
medium, where partons rescatter often, /x 2 is determined by the frequency of 
rescatterings (see Section 3.2). 
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Figure 1: Graphical representation of the 
QCD transport equations (1,2) defining the 
parton cascade model. The 2 — > 2 scatter- 
/i 2 ing process at momentum scale Q 2 is fol- 

lowed by the virtuality evolution from Q 2 
to fl 2 . 

In the leading- logarithmic approximation (LLA), the differential cross sec- 
tion described by the Feynman diagram of Fig. 1 factorizes into a product of 
terms for each of the two incoming and outgoing branch processes and one for 
the binary scattering process. Each branching term, in turn, is represented by 
a product of factors describing the individual branching events and the prob- 
abilities for the partons not to branch further in between. In other words, the 
Feynman diagram of the type shown in Fig. 1 defines a Markov process, as 
it would be expected for any process that can be described by a probabilistic 

"Strictly speaking, Si describes the differential branching probability for an infinitesimal 
change in virtuality; the diagram is an integral representation of Si. 
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one-body transport equation. This statement is no longer valid, if one tries 
to go beyond the LLA. However, certain effects beyond the LLA can still be 
described in terms of conditional probabilities, such as angular or fc^-ordering 
of limited gluons and soft-gluon interference effects. These effects are quanti- 
tatively important andiave been incorporated in parton cascade codes u. 

Xiong and Shuryaklla have argued that one should include all 2 — > n parton 
tree diagrams in the collision term of (|l]), instead of relying on the LLA. Of 
course, if one does this, it is necessary to set the right-hand side of (Q) equal to 
zero, lest one double-counts parton splittings. However, it is important to note 
that parton splittings are related by unitarity to loop diagrams that describe 
the running of the strong coupling constant a s (Q 2 ): 

Both splittings and a s -running are described consistently in the LLA, 
which therefore satisfies the unitarity condition. In plain terms, the combined 
probability for all 2 — > n parton diagrams with n > 3 reduces the probability 
for the occurrence of a 2 — > 2 scattering and so on. By summing all 2 — > n free 
diagrams, but not including the associated loop diagrams, unitarity is violated. 
This leads to the divergence of the sum over n for rather small cm. energies, 
even in the presence of an infrared cut-off for the internal propagators. If 
one naiv&U includes the full sum at energies below this artificial divergence, 
one finds t3 that the chemical equilibration of parton distributions is driven by 
processes with many final-state gluons, i.e. by 2 — > n scatterings with n > 3. 
In my opinion, this result indicates that this approach is questionable. If one 
includes scattering processes with more than two quanta in the final state, one 
also needs to account for unitarity corrections due to loop diagrams. This is 
especially important in the presence of a medium, where loop diagrams are 
enhanced by in-medium contributions (see Section 3.1). 

In the following I will discuss two important issues: 

1. The space-time picture governing the initial-state parton distributions. 
The issue is closely connected with the problem of the decoherence of the 
initial parton wavefunctions. 

2. The problems of infrared divergences of the perturbative parton cross 
sections. This issue is closely related to in-medium corrections to these 
cross sections, as well as to coherence properties of the initial-state wave- 
functions. 

The approach to local thermal equilibrium has heexL extensively studied within 
the framework of the parton cascade picture. Without repeating the 

detailed arguments here, the PCM approach predicts a very short kinetic equk 
libration time, Tth 1 fm/c, which is confirmed by full numerical calculationsH 
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2 Initial state space time picture 

The probabilistic interpretation of parton distributions measured in deep- 
inelastic scattering is based on a summation over all final hadronic states. 
A similar interpretation of one-body distributions arising in transport theory 
is based on the low-level truncation of the (BBGKY) hierarchy of Green func- 
tions and on an expansion in power of h. The validity of this picture ultimately 
relies on the separation of time amies in a dynamic process. Although these 
issues are generally well known E3, their full implications for nuclear parton 
cascades have not been explored. Recent advances Ej in our understanding 
of multiple scattering in QCD have shed some light on the intricacies of the 
formation time concept in non-abelian gauge theories, but more needs to be 
understood. 

The original parton cascade model relied-pn some basic assumptions about 
initial parton distributions in space-time. 113 Denoting the parton light-cone 
momentum by p + , the parton distributions were assumed to be distributed lon- 
gitudinally according to the uncertainty relation: Ap + Ax~ > h. Soft partons 
have the widest distributions in the variable x~ , and were assumed to travel 
both ahead and behind the Lorentz contracted valence quark distributions. 
The argument was that this would not violate causality, because soft partons 
are emitted at a long distance before the collision and, travelling at the speed 
of light, can arrive significantly ahead of the quarks that emitted them. 

The space-time picture of soft partons has been put on a much firmer foun- 
dation in recent years by the work_of McLerran, Venugopalan and others on the 
random light-cone source modelEl (RLSM). In this model, one views the va- 
lence quarks constituting the fast-moving nucleus as a thin, Lorentz contracted 
sheet of locally random color sources. The color source is locally random, be- 
cause valence quarks from several nucleons contribute at the same point in 
transverse space. The area density of color sources is given by /i 2 = 3A/irR 2 , 
where A is the nucleon number and R is the nuclear radius. Clearly fi w A 1 / 6 , 
hence fi becomes a large scale[] for sufficiently heavy nuclei, and a s (fj, 2 ) -C 1 
can serve as basis for a new type of perturbative expansion. Formally, the 
model maps into the problem of weakly coupled QCD in the presence of a 
random two-dimensional color source. 

As shown by KovchegovcJ, this model can be rigorously derived by stan- 
dard light-cone techniques, which permit an explicit representation of the 
Gaussian ensemble of color sources. This representation can also be used to 
calculate the perturbative emission of soft gluons in collisions between two 
nuclei, described as collision between two sheet-like clouds of valence quarks. 

b In practice, < 1 GeV even for the heaviest nuclei. 
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At leading order this soft gluon radiation is given by: 

dN g Aal Nj-\ f 2 F(qa)F(\k-g\a) 

-{1ab{o)) / a g- 



dyd 2 k ± d 2 b n 2 kl N c v At! ^ " J y q 2 (k-q) 2 

where F(qa) is the color dipole form factor of±hc nucleon and TAB(b) denotes 
the nuclear profile function. It can be shown c3 that this classical gluon radi- 
ation smoothly matches onto the perturbative minijet production of gluons at 
higher k±. 

Going beyond the classical approximation by including gluon loop dia- 
grams leads to a better and more rigorous understanding of the space-time 
distribution of soft gluons in a heavy nucleus. c3 The quantum corrections can 
be formulated in the framework of a space-time analogue of the renormaliza- 
tion group equations, describing the cascade of gluon emission leading, to a 
power-law enhancement of soft gluons similar to the BFKL equation. E3 The 
RLSM approach also describes saturation effects in the parton distribution at 
small fcx • It would be very .interesting to study its relation to the perturbative 
theory of gluon shadowing in this regime. 

The picture that emerges is the following: Gluons in the classical field 
generated by the valence quarks are fully Lorentz contracted by the Lorentz 
factor 7 associated with the colliding nuclei, but gluons spawned by a splitting 
of those primary gluons experience only a partial Lorentz contraction of order 
ary, where x is the momentum fraction carried by the parent gluon. As the 
branching process evolves to softer and softer gluons, the spatial extent of this 
gluon cloud becomes more and more diffuse in the light cone variable x~ . This 
result confirms the intuitive picture embodied in the original parton cascade 
model, and provides a quantitative formulation of it. 



3 Medium effects 

Medium-effects on parton-parton interactions are essential for the viability of 
the parton cascade model. To wit, the application of the PCM to nucleon- 
nucleon collisions requires the introduction of ad hoc cut-offs describing non- 
perturbative QCD effects, such as quark confinement and chiral symmetry 
breaking. Medium effects, which grow rapidly in size as function of A, can 
produce perturbative cut-offs when the density of the medium becomes suffi- 
ciently high. E.g., QCD is known to become perturbative [] at high temperature 
when the electric screening mass fin 3> Aqcd- 

c Some nonperturbative effects remain even at high T, precisely because static magnetic 
interactions are not screened by perturbative in-medium interactions. 
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Thus, medium effects work in favor of the parton cascade model. The 
problem is that medium effects are very complicated and not easily treated 
correctly. The two main medium effects that are known to provide effective 
infrared cut-offs to perturbatively divergent parton interactions are: 

• (color-electric) screening, which suppresses soft 2 — > 2 scattering ampli- 
tudes; 

• (gluon) radiation suppression, which dampens 2 — > 3 (and 2 — > n) am- 
plitudes with soft particles in the final state. 

Dynamical screening, at lowest order, is described by the in-medium contri- 
butions to the one-loop gluon polarization function: 




At moderately high q±, the gluon population grows like n(k) w (A1A2H/ 3 , 
providing a screening scale n{A) that increases rapidly with nuclear size. E3 

Radiation suppression, also known as the Landau-Pomeranchuk-Migdal 
(LPM) effect, is a much more complicated effect. Its theoretical description 
requires a good understanding of the multiple scattering problem in QCD. 
Considerable progress has-cecently been_made in this_area, especially through 
the work of Baier et al. c2, ZakharovEil and others E3. The main difference 
between QCD and the well-known case of QED is that a radiated gluon also 
rescatters in the medium at the same order in a s as the emitting particle; this 
is not so for a photon in an electromagnetic plasma. 

Revisiting the diagram shown in Fig. 1, the parton cascade model requires 
infrared cut-offs for both the central 2 — > 2 scattering matrix element and each 
of the four branching cascades. This is where medium effects help: at high 
density one expects the in-medium modifications of the elementary scattering 
amplitude to ensure an infrared safe behavior. To date, two attempts have 
been made to practically implement the action of these medium effects: 

(1) In the self-screened parton cascade modelEl (SSPCM), the color-electric 
screening scale ^{pt) was calculated self-consistently for primary parton 
interactions only, and the further evolution of the parton plasma was 
described in the framework of the hydrodynamical model. 

(2) In the PCM code VNlS, the space-time picture of parton interactions 
is linked with the virtuality evolution of partons. A new interaction is 
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permitted if its momentum transfer exceeds the virtuality scale of the 
participating partons at that time. Soft-radiation is suppressed in the 
medium by this rule, because intermediate collisions continuously reset 
the parton virtuality. 

Both of these approaches are based on specific assumptions about the relation 
between the space-time and virtuality evolution of off-shell components of the 
parton distributions. One way to study this issue vigorously is the Wigner 
function representation. First results EJ obtained by this method are inter- 
esting, but do not yet fully address the complications encountered in a QCD 
parton cascade. A different approach to the problem of the space-time picture 
of off-shell quantum fluctuation is based on a modification of the QCD evolu- 
tion equationsEj to include an infrared scale. In free space this infrared scale is 
determined by properties of the final state (hadronization scale) ; in a medium 
it is determined by screening effects. 



3. 1 Self-screened parton cascade 

Here one considers the scattering of an initial state parton as completed after 
a time t(pt) which depends on the momentum transfer in the reaction. The 
uncertainty relation suggests t(j?t) ~ h/pT- (We will drop the factor ti in the 
following.) The scattered partons are then assumed to screen the scattering 
processes that involve a smaller momentum transfer: 

3 f°° 
Md(Pt) = — a s {p 2 T ) / d 3 k\\7 k n(k)\. 

The density of partons scattered at pr is, in turn, influenced by the screening 
because the differential cross section depends on /j,q' 

da_ a s {p T ) 2 | , 1 2 

d P 2 T ~ (p 2 t + vUptW 1 [ ' j| ■ 

If ht>{jpt) becomes large enough at low p?, so that da j dp\ remains perturba- 
tively small, the coupled set of equations can be integrated down to px = 0. 
Since the rapidity density of scattered partons grows as (A1A2) 1 / 3 ln-^s, this 

sely 
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condition requires large A and high energy. The SSPCM concept is cL 
related to the RLSM approach proposed by McLerran and Venugopalan : 

Quantitatively, one finds that /id approaches about 1 GeV at low px in 
Au + Au collisions at RHIC energy (100 GeV/u) and 1.5 GeV at LHC energy 
(2.75 TeV/u). The differential minijet cross section as function of pr peaks 
at about the same value, clearly showing the improved infrared behavior of 



S 



the self-screened parton cascade. The total deposited energy within one unit 
of rapidity and after a characteristic formation time of 0.25 fm/c is eo ~ 60 
GeV/fm 3 (KHIC) and e w 430 GeV/fm 3 (LHC). The conditions established 
by the SSPCM can, therefore, be taken as initial conditions for the thermal 
and chemical evolution of a quasi-equilibrated parton plasma. 

The equations for the evolution p£ such a plasma were formulated by Biro 
et al. c3 and by Xiong and ShuryakOJ. Extensive calculationsBEa, including 
longitudinal and transverse expansion, have shown that the plasma cools down 
to the critical temperature of QCD (T c ta 150 MeV) after 5 fm/c (RHIC) and 
10 fm/c (LHC). The emission of electromagnetic probes by such an evolving 
QCD plasma has also been calculated. I?J 

3.2 Monte-Carlo space-time cascade 

The statistical implementation of a parton cascade by means of a Monte-Carlo 
code (VNlEa) achieves an improved infrared behavior through heuristic rules 
that suppress certain interactions on the basis of kinematic considerations. 
The first rule asserts that independent scattering events involving the same 
parton require a sufficient time separation so that the time between scatterings 
is larger than the duration of the individual events. With the duration of an 
interaction again defined as t(pt) ~ Vt 1 where pt is the momentum exchange, 
this requires that the time between interactions At > t{j>t). Another way of 
ensuring this condition is to endow a parton after a scattering by pt with 
an initial virtuality Qo — pt, which then gradually decreases with time as 
Q{t) = Qqt(pt)/t. A subsequent scattering with p' T requires that p' T > 
Q(t) at the moment of the interaction. A second similar rule suppresses soft 
parton splittings in the presence of multiple scattering. Again, this rule can 
be formulated in terms of a parton virtuality that decreases with time between 
scatterings and is reset by each new interaction (see Fig. 2). 

Figure 2: Illustration of in-medium sup- 
pression effects incorporated in the VNI 
code. The virtuality of a scattered parton 
evolves with time, Q(t). Sequential scatter- 
ings or branchings are suppressed, if Q(t) 
is too large. 

Although VNI still contains "arbitrary" infrared cut-off parameters (deter- 
mined by a comparison with nucleon-nucleon interactions), these are needed to 
limit soft scatterings in the initial set of parton interactions, and they become 
important toward the end of the cascade evolution when the parton plasma 
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becomes more and more dilute. The cut-off parameters effectively determine 
the end of the cascade evolution and, by suppressing soft interactions, they 
trigger the hadronization of parton clusters. 

The fact that this improved parton cascade is much less dependent on arbi- 
trary cut-off parameters makes it possible to explore applications of this model 
at lower energies where nucleonpnucleon collisions provide little information. 
Recently, Geiger and Srivastavac9 have studied the predictions that the VNI 
code makes for nuclear collisions in the energy regime of the CERN-SPS. While 
most of the particle yield at rapidities \y\ > 2 is produced by fragmentation of 
the unscattered beam remnants, the model predicts a significant contribution 
to particle production at r central rapidity from partons that have undergone 
perturbative scattering. EH This contribution is rising rapidly with nuclear 
mass A, roughly as (A1A2) 1 / 3 . This perturbative contribution to the energy 
deposition at \y\ < 1 coincides with a rapid increase of the energy density in 
scattered partons at r < 1 fm/c, which rises from about 2 GeV/fm 3 in S+S to 
5 GeV/fm 3 in Pb + Pb. This rise may be correlated with the much enhanced 
suppression of charmopium production in Pb + Pb collisions as observed by 
the NA50 experiment. e3 

4 Summary 

The parton cascade model was developed to provide a QCD-based description 
of the approach to a locally thermalized state in collisions of heavy ions in the 
RHIC energy regime and beyond. In its original formulation the PCM predic- 
tions were critically dependent on several cut-off parameters that had to be 
determined from pp collision data. Recent advances in incorporating medium 
effects into the parton interactions have reduced this dependence significantly, 
possibly allowing the application of the PCM over a wider energy range. Re- 
sults obtained for nuclear collisions at CERN-SPS energies are intriguing. 

The in-medium effects that modify parton-parton interactions not only 
reduce the parameter dependence of the model, they also provide valuable 
insight into the dynamics of a dense parton plasma. It is clear that we are 
here just at the beginning. The transport properties of off-shell quanta need 
to be understood much better, not only in cases where the off-shell propagator 
is dominated by a well-defined resonance, but especially in the case where the 
particles never get close to their mass shell as it applies to QCD. Another 
open question concerns the need for mean color fields. Such fields are not 
included in present versions of the PCM, but the random light-cone source 
model suggest that mean fields may be essential ingredients of a complete 
description of soft processes in nuclear collisions. One may have to take an 
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average over a Gaussian ensemble of mean fields, thus the width of the field 
distribution may be more important than the expectation value which remains 
zero. It would be interesting to explore possible_cprmections of the RLSM to 
the traditional chromo-hydrodynamical model. c3o 

Ultimately, the question is whether the parton cascade model can be re- 
placed by a controlled approximation scheme where, the principle, successive 
orders of ever more sophisticated corrections can be calculated. We are still far 
away from a consistent formulation of transport phenomena off equilibrium in 
QCD, and it is even unclear whether we have a full insight into what are the 
small parameters in such an approximation scheme. It is clear that a high den- 
sity of excitations of the QCD vacuum is an essential condition, but there are 
many subtleties if one wants to go beyond this statement. However, progress 
in this field is steadily made, and there is reason to hope that a consistent 
formulation of transport phenomena in QCD can ultimately be achieved. 

My final remarks concern the treatment of the late phase of a heavy ion 
collision when the dense matter breaks up into individual hadrons. Here one 
has the choice of either applying a phenomenological hadronization model to 
make the transition to a hadronic cascade, or to change from the parton cascade 
model to a hydrodynamical description already in the dense plasma phase as 
soon as approximate local (kinetic) equilibrium has been reached. The second 
scheme has the advantage that an equation of state describing the QCD phase 
transition can be easily incorporated; it has the disadvantage that another 
transition to some cascade-like scheme is needed in the very late phase when 
the mean free paths of particles become too long to sustain the hydrodynamical 
description. 
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